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Abstract 

The time independent spherically symmetric solutions of General Relativity (GR) coupled to a 
dynamical unit timelike vector are studied. We find there is a three-parameter family of solutions 
with this symmetry. Imposing asymptotic flatness restricts to two parameters, and requiring that 
the aether be aligned with the timelike Killing field further restricts to one parameter, the total 
mass. These "static aether" solutions are given analytically up to solution of a transcendental 
equation. The positive mass solutions have spatial geometry with a minimal area 2-sphere, inside 
of which the area diverges at a curvature singularity occurring at an extremal Killing horizon that 
lies at a finite affine parameter along a radial null geodesic. Regular perfect fluid star solutions are 
shown to exist with static aether exteriors, and the range of stability for constant density stars is 
identified. 
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I. INTRODUCTION 



For nearly a century Lorentz invariance has been a basic assumption in theoretical physics. 
However, various approaches to the problem of quantum gravity question whether Lorentz 
symmetry is truly fundamental. In this context, it is useful to consider low energy Lorentz 
violating (LV) models from the point of view of effective field theory. It makes sense to 
distinguish LV effects in the matter sector from LV effects in the much more weakly coupled 
(at low energies) gravitational sector. LV effects in the matter sector are tightly constrained 
by h,gh precision tests of LorenU invariance fl. Lorentz vio.afon in the gravitational 
sector was studied by Will and Nordvedt [2|] in the early 1970's, and it was pursued in 
the 1980's by Gasperini [s] motivated by singularity prevention, and by Kostelecky and 



Samuel 



J] motivated by the possibility of spontaneous LV in string theory. More recently a 



number of approaches to incorporating LV into gravity have been explored; see for example 



Refs. 



ler oi approac 
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ll| and references therein. 



Here we study the theory of a unit timelike vector field coupled only to gravity, as an 
LV modification of General Relativity (GR). To preserve general covariance this field must 
be dynamical. It can be thought of as the 4- velocity of a preferred frame, which breaks boost 
symmetry at every point of spacetime. This property it has in common with the 19th century 
concept of the aether, and like that aether it is also a dynamical system. Hence a suitable 
name for it is the "aether" , although it has nothing to do with electromagnetism. Since the 
aether is coupled to Einstein GR we refer to the theory as "Einstein-Aether" theory; we 
sometimes use the nickname "ae-theory" for short. For a review of the motivation, history 
and recent status of the Einstein-Aether theory, see |l^ and the references therein. 

Observational constraints on the Einstein- Aether theory have been determined from PPN 




19|, primordial nucleosyn- 



15|.stabi.,tyandene™ypos,t,vityQ.Qfl 
thesis [20|, and Cerenkov radiation 5l8l|. The combined constraints from all of these are 
reviewed in Ref. 15|, and constraints from radiation damping are discussed in 2l|]. Also 
important for determining theoretical viability and observational constraints are the prop- 
erties the spherically symmetric static solutions. These solutions were previously obtainec 
in the the special case where the aether dynamics is Maxwell-like {4, 6|. It was shown in 6] 
that the Reissner- Nordstrom metric in a spherically symmetric static gauge with fixed norm 
is a solution, and it was claimed (incorrectly, as shown here) that this is the only solution 
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in that special case. The asymptotic weak field limit of the general case was studied in 



Ref. 



13| . where it was found that there is a two-parameter family of asymptotically flat 



spherical, static solutions. A thorough examination of the fully nonlinear solutions has not 
been carried out before. 

In this paper and its companion (devoted to black holes) 22(] we complete a general 
survey of the time-independent spherically symmetric solutions. In section [TTl of this paper 
the action and field equations for the Einstein-Aether theory are reviewed. Specializing to 
stationary, spherical symmetry, section [Till studies the structure of the field equations as 
ordinary differential equations and shows that locally there is a three parameter family of 
vacuum solutions. Imposing asymptotic fiatness reduces the number of parameters to two. 

In section HVl the general solution in which the aether vector is aligned with the timelike 
Killing field is found. Being a unit vector, the aether is completely determined by the 
metric in this case. This solution is asymptotically fiat and described by one free total mass 
parameter. For negative total mass there is a naked singularity at the origin. The positive 
mass solution has a wormhole-like spatial geometry, reaching a minimum area 2-sphere at 
some radius like the Schwarzschild solution. Unlike in the Schwarzschild solution this throat 
is not on a horizon. Inside the throat the spheres re-expand to infinitely large size in finite 
affine parameter along a radial null geodesic and finite or infinite proper distance depending 
on the coupling parameters in the Lagrangian. When the distance is finite the internal 
infinity is singular, and it occurs at a would-be extremal Killing horizon. When the distance 
is infinite the metric is asymptotically singular. 

In section |V] it is first shown that pure aether stars do not exist, i.e. there are no 
asymptotically fiat self-gravitating aether solitons with a regular origin. (Given the results 
of the previous section, regular aether stars could only possibly have existed if the aether 
were not everywhere aligned with the Killing vector.) Next it is shown that in the presence 
of a perfect fiuid, regular asymptotically flat star solutions exist and are parameterized (for 
a given equation of state) by the central pressure. For the case of constant density the star 
solutions are found by matching numerical integration for the interior to the vacuum solution 
found in section IIVI As in GR, for a given density there is a maximum mass. Utilizing the 
critical behavior of the mass as a function of stellar radius R it is shown that if they are 
stable at small mass, these stars are unstable beyond the maximum mass. We conclude in 
section |Vl] with a discussion of open questions raised by these results. 
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II. EINSTEIN- AETHER ACTION 



The action for Einstein-Aether theory is the most general diffeomorphism invariant func- 
tional of the spacetime metric Qab and aether field u"' involving no more than two derivatives, 

S = [ L d^x (1) 



where 

L = -R- K'^'mnVau'^Vtu'' - X{gabU^u' - 1). (2) 
Here R is the Ricci scalar, K"'^rnn is defined as 

i^"'™n = cig^'^'g^n + 02^5^ + 035^(5^ + c^u'^u^'g^n (3) 

where the q are dimensionless constants, and A is a Lagrange multiplier enforcing the unit 
timelike constraint. This constraint restricts variations of the aether to be spacelike, hence 
ghosts need not arise. A term of the form RabU°'u^ is not explicitly included as it is propor- 
tional to the difference of the C2 and C3 terms in ([1]) via integration by parts. The metric 

signature is (H ) and the units are chosen so that the speed of light defined by the 

metric gab is unity. In spherical symmetry the aether is hypersurface orthogonal, hence it 
has vanishing twist Ua = eabcdu'^^'^u'^- When u'^ is a unit vector the square of the twist is a 
combination of the Ci, C3, and C4 terms in the action ([1]), 

The C4 term can thus be absorbed by making the replacements ci ci + C4 and C3 — > C3 — C4. 

Following the observational constraints we assume here when studying the fluid star so- 
lutions that the only significant coupling of u°' to matter is through a universal "matter 
metric" g^^^^"^^ = gab + crUaUb, where a is a constant. Replacing gab by g^b^^*'^^ as the in- 
dependent metric field in the action returns an action with the same form as but with 
new values of the constants Ci^2,3,4 that depend on a |23j. Hereafter we assume that such 
a field redefinition has already been performed, so that gab is the metric to which matter 
couples universally. The absence of any other coupling of m° to matter has no theoretical 
justification in this purely phenomenological approach, and may be regarded as unnatural. 
However our goal here is just to explore consequences of gravitational Lorentz violation in 
a phenomenologically viable setting. It remains an open question whether this can emerge 
as an approximation to a more fundamental underlying theory. 



The field equations from varying ([T]) plus a matter action (coupled only to the metric) 
with respect to g"'^, u"" and A are given by 

Gab = T^^ + SvrGT^,, (5) 

Va - CiilaSmU" = XUm, (6) 

gabU'^u' = 1, (7) 

where 

= K'^'mnVbU'^. (8) 

The aether stress tensor is given by 13| 



+ci [(V^Ma)(V™Mb) - (VaM™)(Vfen")] 

+C4 ilallb 

—Lugab, (9) 

where L„ = — i^"''mnVaM™'VbM". The Lagrange multiplier A has been eliminated from ([9]) 
by solving for it via the contraction of the aether field equation (j6]) with u"". 

Some words about terminology are in order. Spacetimes admitting a timelike Killing 
vector field are generally called stationary. In the special case where is hypersurface 
orthogonal, and therefore invariant under a time reflection t —t, the spacetime is said to be 
static. A stationary aether field on a stationary spacetime is one whose Lie derivative with 
respect to vanishes. If the spacetime is static, one might be tempted to say the aether is 
"static" , however this is not really appropriate since the aether itself breaks the Killing time 
reflection symmetry. The solutions studied in this paper involve a static metric coupled to a 
stationary aether. This general situation will be called here "stationary spherical symmetry" . 
An important special case occurs when the aether is parallel to the Killing vector. We refer 
to this special C3iSG clS "static aether". Such an aether changes sign under the Killing time 
reflection, however the action ([T]) is invariant under —>■ —u"" so the sign of u"" has no 
physical meaning. Note that regular black holes cannot have static aether fields since the 
Killing vector is null, not timelike on the horizon. 
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III. CLASSIFICATION OF STATIONARY SPHERICAL SOLUTIONS 



Stationary spherically symmetric solutions describe, for example, a black hole or the exte 



24| 



rior of a time-independent star. In spherical symmetry all stationary metrics are static 
The line element can be written in Schwarzschild type coordinates, 

ds^ = e^'-'-Ut^ - B{r)dr^ - r^dQ^, (10) 

and the aether field takes the form 

u = a{r)dt + h{r)dr. (11) 

The unit constraint on becomes 



e 



^Wa(r)2 - B{r)h{rf = 1, (12) 



which can be used for example to eliminate b{r). The t-component of the aether field 
equation ([6]) can be used to solve for A in this case, and the remaining field equations 
reduce to five ODE's: the tt, rr, tr, and 66 components of the metric field equation and the 
r component of the aether field equation. These five equations involve the eight functions 
{A", A', A, B', B, a", a', a}, where prime denotes differentiation with respect to the argument 
r which is suppressed. The equations are too complicated to be worth writing down here, 
so we shall just describe their structure. Using the tt and 66 metric equations along with 
the r component of the aether equation, one can solve for A", a", and B' in terms of the 
remaining five functions {A',A,B,a',a}. It turns out that only one additional piece of 
information remains in the tr and rr equations, which can be used to solve (for example) 
for B in terms of {A', A, a', a}. Finally, v4(ro) at any given value r = rg can be chosen 
at will by allowing for an appropriate scaling of the t coordinate. At a given rg value, the 
remaining three values {A'(ro), a'{rQ), ci(ro)} then determine a (local) solution by integration 
with respect to r. This shows that there is in general a three-parameter family of spherically 
symmetric stationary solutions. 

To illustrate this reasoning in a more familiar setting, we can apply it to the field equations 
of pure GR in Schwarzschild coordinates, 

Gtt oc rB' -B + B^ = (13) 
Grr oc rA' -B + 1= (14) 
Gee oc 2rA"B + A'{2B - rB') + rA'^B - 2B' = 0. (15) 



These can be used to solve for B\ B, and A" in terms of A' and A. Using the freedom to 
scale t the initial value A{ro) can be fixed at will, so we recover the well-known fact that 
static spherically symmetric solutions in GR are characterized by one free parameter, in this 
case the value of A'{rQ). In the Einstein- Aether theory the aether vector and its derivative 
provide two additional degrees of freedom at each point. 

Birkhoff's theorem in GR states that the only spherically symmetric solution is static 
and given (up to coordinate freedom) by the Schwarzschild metric. The radial tilt of the 
aether provides another local degree of freedom in ae-theory, so spherical solutions need 
not be time-independent. But, as we have seen, even restricting to stationary spherically 
symmetry ae-theory has more solutions. In this paper we will focus primarily on the static 
aether solutions, which form a one-parameter family. Black hole solutions, which comprise 
a different family, are studied in a companion paper. 

Numerical integration of the ae-theory field equations as ODE's out from some arbitrary 
point ro with generic initial conditions yields singularities in A{r), B{r), and a{r). However, 
there is a two-parameter family of asymptotically flat solutions. This was first found in 



Ref. [13|] using a perturbative expansion about infinity. Asymptotic flatness was imposed 
there by assuming regular power series expansions about x = 1/ p = 0, where p is the 
isotropic radial coordinate. Asymptotic flatness can also be imposed using the "shooting 
method". This is simple to implement here since it is only necessary to tune one of the 
three initial values {y4'(ro), a'(ro), a(ro)} so that, for example, A{r) approaches a constant 
value as r ^ oo. The field equations then automatically enforce the remaining asymptotic 
flatness conditions. In GR, by contrast, asymptotic flatness is a consequence of the vacuum 
fleld equations without any tuning of initial data, so the one-parameter family of local 
(Schwarzschild) solutions is automatically asymptotically flat. 

IV. STATIC AETHER 

In this section we obtain the static aether solution, where the aether vector is propor- 
tional to the timelike Killing fleld C,"" and therefore entirely determined by the metric. 
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A. Field equations with static aether 

Using the Schwarzschild type coordinates in (fTO!) and (fTT]) . the static aether has b{r) = 
and a(r) = exp(— y4(r)/2), i.e. 

u = (16) 

In this case the C2 and C3 terms drop out of the field equations. To see why, note that 
( IT6|) imphes VaU"" = 0, so all variations of C2 term in the action ([1]) vanish. In addition, 
the normalization UaU"' = 1 implies UbVaU^ = 0. These conditions together with spherical 
symmetry imply that the derivative of u'' has the form 

"^aU' = UaS\ (17) 

where is a radial vector orthogonal to (We note in passing that contraction of (fT7|) 
with u°' reveals that is the acceleration of the aether worldlines.) Therefore (VaM^)(VfeM") 
vanishes, so the variation of the volume element in the C3 term of the action ([T]) vanishes. 
The remaining variation of the C3 term is proportional to UaS^diVhu"') = s^6{uaVbU°'), which 
vanishes for all variations {6gab,Su°') preserving the normalization gabU°'u^ = 1- Moreover, 
as explained in Section [Tll in spherical symmetry the C4 term can be absorbed into the Ci 
and C3 terms, hence the solutions with static aether are fully characterized by the case with 
only Ci non-zero. 

When only Ci is nonzero the aether field equation (P) reduces to 

CiVVaM^ = Xu\ (18) 

Using ( ITTI) this becomes 

CiM'^VaS* = Xu\ (19) 

Contraction of the left hand side of f|T9|) with Sb is proportional to w^VaS^, which vanishes 
since is a scalar that must be constant along the Killing direction parallel to m". Therefore 
both sides are parallel to m*, so the aether equation only determines A. Contracting both 
sides of (IT^ with Ub we find 

A = -cis^ (20) 
having made use of (ITTl) and u^VaSb = —s^VaUb, which follows from u'^Sb = 0. 
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The metric field equation is Eab = Gab — T^b — 0; '^^i components of 

Eab are given by 

Ett = (e^/r'^B) + B + rB'/B) - iy{8rA' + r^A'"^ - 2r'^A'B'/B + Ar^A")] (21) 
Err = r-2(l - 5 + rA' + ur^A'^) (22) 



Eaa = B ^ 



{2rA' - 2rB'/B + r^A'^ - r^A'B'/B + 2r^A")/A - z/rM'^ 



(23) 



wliere for notational convenience we liave introduced tlie symbol 

Cl 
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(24) 



Using the Err equation one can solve for B, 



B = 1 + rA' + ur^A''^ 



(25) 



Substituting this solution for B into the Eu and Err equations, we find that the equations 
are redundant and the system is described by the second order ODE 



r^A" + 2rA' + r'A''' + v r" A' = 



1' , „2 A/2 



3 Af3 



(26) 



A constant shift of A can be absorbed by a scaling of the t coordinate, hence there is 
just a one parameter family of solutions. As in GR, the solutions in this family are all 
asymptotically flat. 

B. Static aether solutions: general analysis 



To solve (|26|) we define the function Y{r) by 



Y = rA\ 



in terms of which the solution for B becomes 



B = 1 + Y + uY' 



(27) 



(28) 



and Eqn. fl26l) for A becomes 



dY/dr = -(F/r)(l + Y + uY^) 



(29) 
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The problem is thus reduced to quadratures: integration of this equation yields Y{r), which 
also directly yields B via fl28l) . To determine A we combine (1271) and fl29l) to obtain 

dA/dY = -1/(1 + Y + iyY^), (30) 

which yields A{Y) by integration. 

The character of the solutions is evidently controlled by the roots of B, 



in terms of which we have 



Kt = (-1 ± VT^)/{2u), (31) 



B = iy{Y -Y_){Y -Y+). (32) 



The nature of the roots depends on the value of u. We consider here only positive u, since 
that is required by positivity of the energy of linearized spin-0 waves |19], and we restrict to 
z/ < 1/4 since the Newton constant = G/{1 — C14/2) becomes negative beyond this limit. 
In the pure Ci case also the stability [igI] or positive energy [l^ of linearized waves requires 
Ci < 1, or < 1/8. One can visualize the roots graphically: they are the intersections of 
the line Y + 1 with the inverted parabola —uY"^. When z/ = 1/4 the parabola is tangent to 
the line, and the two roots coincide at Y = —2. The larger root approaches — 1 as z/ ^ 0, 
while the smaller root approaches —00, hence in the range < z/ < 1/4 the roots fall within 
the ranges 

- 00 < F_ < -2, -2<Y+< -1. (33) 

Note that Y_ = l/(z/F+), and z/ = -(1 + F+)/F^. 

We can integrate ( 130|) and ( 129|) to find both A and Y using the factorization ( l32i) and 
partial fractions. The result is 

F/F_ 



and 



1 



Y \ fY -Y\^ . . 

where rmin is an integration constant. The graph of r/rmin vs. Y is plotted in Fig. [H for the 
case Ci = 1. The values of the sphere radius r and metric functions B and at the special 
values of Y are given in Table [H 
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FIG. 1: Graph of r/rmin vs. Y for ci = 1. The curves approach 1 asymptotically on both sides. The 
range (1 + , 0) defines a negative mass solution with naked singularity at y = 1+ and asymptotically 
flat region at y ^ 0. The range (0, oo) defines a positive mass solution, with a minimal 2-sphere as 
Y — > oo. The range (— oo, YS) continues that solution to the other side of the minimal sphere, with 
a singularity at a sphere of infinite radius at Y = Y-. This sphere lies at finite radial distance if 
ci < 3/2. There is no solution with timelike aether in the range (y_, y+) since the radial coordinate 
is timelike there (see section HV B 4p . 

TABLE I: Sphere radius and metric functions at special Y values. 



Y 


y_ 


y+ 





iboo 


r 


oo 





oo 


'"min 


B 








1 


oo 


N 





oo 


1 


> 



1. The GR limit: Schwarzs child solution 

To help to interpret the general case, we consider first the pure GR limit ci = 0, for which 
y+ = — 1 and y_ = — oo. The solution is then 

B = l + Y (36) 
N = 1/(1 + Y) (37) 
Wr = Y/{l + Y). (38) 

This is just the Schwarzschild solution, with Y = l/{r — r^i^) and rmin = 2M. Spatial 
infinity corresponds to y = 0, and as y ^ oo the radius decreases to rmin at the bifurcation 
surface of the horizon. The other side of the wormhole is here labelled by the same values 
of y. The range — oo < y < — 1 corresponds to the future wedge of the black hole interior, 
where the Killing vector is spacelike. The remaining range — 1 < y < is also significant. 
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It corresponds to the negative mass Schwarzschild solution. 



2. Static aether solutions for generic ci 

For generic values of ci the limit F — still corresponds to an asymptotically flat spatial 
infinity, where the limiting form of the solution is 

B = ! + ¥ + ■■■ (39) 

N = !-¥ + ■■■ (40) 

Y = 2M/r + ■ ■ ■ , (41) 

and the mass M is related to the minimum radius by 

W2M = i-Y^)-\-l - (42) 

This ratio grows smoothly from 1 for Ci = 0, to about 1.23 for Ci = 1, and reaches e/2 ^ 1.4 
for Ci = 2. 

Series solution in powers of a; = 2M/ r yields 

B = l + x + (l + u)x^ + ■■■ (43) 

N = 1-x- (z//6)x^ + • • • (44) 

Y = x + x^ + {l + u/2)x^ + ■■■ . (45) 

To leading order in 1/r this agrees with the Schwarzschild solution, as already seen previously 



m [13| , where the two Eddington-Robertson-Schiff PPN parameters 7 and (3 were found to 
agree with the GR value of 1. This completes our characterization of the asymptotically fiat 
region. What happens when we follow the solution to smaller values of r? 

The answer depends on the range of Y considered. For Y G {Y+,0), equation fl35|) or 
its graph in Fig. [1] indicate that rmin niust be negative, which according to (H2l) implies a 
negative total mass M < 0. In this case there is a naked singularity at Y = Y^ {r = 0) 
connected to an asymptotically fiat region at y = (r = 00), like in the negative mass 
Schwarzschild solution. 

For positive Y the solution is different. It is seen again from fl29l) or its graph that 
Y grows monotonically as r decreases. As F — 00 the r.h.s. of (l35l) goes to 1, so this 
limit for Y corresponds to a minimum radius fmin, just as in the case of the Schwarzschild 
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solution. However, the solution behaves quite differently from Schwarzschild. First, instead 
of A^(rmin) = we have 

Y \ -Y+/(2+Y+) 



iV(r„.in) = i^^j , (46) 

so the minimal 2-sphere does not sit at a Killing horizon. The value of A^(rinin) grows 
smoothly from for Ci = to about 0.083 for Ci = 1 and reaches ~ 0.135 for Ci = 2. 
(Recall that in the GR limit we have Y_ — oo.) Another difference due to the finiteness 
of YL is that the solution continues with negative Y values, with the two values Y = ±oo 
identified. According to ( |29ll . as Y grows from —oo up to YL, r increases from rmin to 
oo. Therefore the "interior" of the minimal 2-sphere flares out to infinite radius as in the 
Schwarzschild solution. But unlike the Schwarzschild case, now two values of Y correspond to 
each r, and the "interior" geometry is not equivalent to the exterior. In fact the difference 
is quite dramatic: at the internal infinity both and B go to zero, whereas they both 
approach one in the asymptotically flat region. 

The Carter-Penrose diagram for this solution is the square diamond in Fig. [2]with asymp- 
totically flat past and future null infinity on the lower and upper edges bounding the right 
hand side and a singularity on both edges S"^ and at 5° bounding the left hand side. S"^ are 
singular Killing horizons with vanishing surface gravity. The proper distance to 5**^ along a 
constant t surface is finite if < Ci < 3/2 and infinite if 3/2 < ci < 2, while for any ci the 
affine parameter to along a radial light ray is finite, as we now demonstrate. 

As the internal infinity at K. is approached, the proper radial distance on a constant t 
surface behaves as 

f - = 5^/2 _ _ y)l/2_ (47) 

dr 

On the other hand (|35|1 shows that in this limit the relation between Y and r is 

(r_-F)~r^, (48) 

so we have 

dl 

— = r^+^y+. (49) 
dr 

The exponent of r in (H9l) is always negative, and it is equal to —1 when y+ = —4/3, which 
corresponds to u = 3/16, i.e. ci = 3/2. For ci < 3/2 the radial distance to is finite and 
there is a curvature singularity at r = oo that shows up, for example, in the square of the 
Riemann tensor. 



13 



FIG. 2: Carter-Penrose diagram of the static aether solution. The left hand edge corresponds to 
spheres of infinite radius and is singular. 

Along a radial null geodesic the quantity Nt = Ndt/dX is conserved, where A is an affine 
parameter. Together with the lightlike condition Nt"^ — Br"^ = this implies that as the 
internal infinity is approached the affine parameter behaves as 

dX/dr = {NBf^ ~ (y_ - r)V(2+y+) ^ (50) 

The affine parameter distance to S"^ is therefore finite for all E (—2, —1), corresponding 
to all ci e (0,2). 

Note that since the minimal 2-sphere is not hidden by a horizon, a spherical congruence 
of null rays will converge towards the minimal sphere and exit the other side with a positive 
expansion. The Raychaudhuri equation shows that this can happen only if Rabk"'k^ < 
somewhere along the congruence, so we infer that the aether stress tensor must violate 
the null energy condition in this solution. We computed the curvature for this solution 
and found that Gu = —uY'^/Br'^, so the energy density of the aether (oc Gu) is negative 
everywhere . The solution nevertheless has positive total mass, which may at first seem to 
be inconsistent but it is not. The total mass of an asymptotically fiat spacetime is given 
by a surface integral at spatial infinity, which for stationary spacetimes is proportional to 
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the volume integral j^Rabn°'^^dV , plus a surface term if there is an inner boundary 24 1. 
Since Rab oc Tab — {^/2)Tgab, it is not the energy density that figures in the total energy 
but rather Ru. Quite surprisingly, it turns out that Ru vanishes everywhere in the static 
aether solution. (The only nonzero component of the Ricci tensor is Rrr-) Hence the 
energy integrand vanishes identically, as in Schwarzschild spacetime. This does not mean 
the total energy vanishes however, since there is a contribution from the inner boundary. In 
Schwarzschild that inner boundary may be pushed off to the asymptotic region on the other 
side of the Einstein-Rosen bridge, but the static aether solution is singular on the other side 
of the throat. One can think of the mass as determined by a boundary condition at this 
singularity. 

Let us briefly consider solutions for ci in the range 3/2 < ci < 2, corresponding to 
3/16 < z/ < 1/4. In this case the distance to the internal infinite radius sphere is infinite, and 
all the algebraic and differential invariants that we checked (including R, RabR'^'^, Rabcd.R"''^'^'^ , 
RabU°'u^, and {V aU^){Sbu"')) are asymptotically zero. However, the curvature component 
Rabk°'k^ blows up asymptotically, when fc" is the tangent to an affinely parameterized radial 
null geodesic approaching the internal infinity. The invariant u"'ka blows up as A^~^/^, since 
j\j-i/2^a jg Killing vector and k°'^a is conserved along the geodesic. This suggests that 
the above-mentioned invariants vanish because the tensor structure of the curvature, the 
aether, and all derivatives is determined by a single null vector pointing in the future radial 
null direction opposite to A;", i.e. pointing away from the internal area-infinity. 

Returning now to the generic solution for < Ci < 1, we examine more closely the 
behavior at the throat and at the internal infinity. Since the spherical radius r is not a good 
coordinate at the minimal area sphere, nor at the internal area-infinity, we adopt instead 
the proper length coordinate /, in terms of which the line element takes the form 

ds^ = N{l)df - df - r{lfdVL^ (51) 

To get an idea of how the throat geometry depends on ci, we fix the mass M of the solution 
and plot in Fig. [3] the numerically computed function r(/) for several different values of Ci. 
There is a discontinuity at ci = where the solution abruptly changes from a singular fiare- 
out in finite proper distance to a perfectly regular Einstein-Rosen bridge. The singularity 
approaches the throat as ci —^0, but in the same limit the curvature becomes finite and the 
other half of the bridge suddenly appears. 
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FIG. 3: Plot of area radius r vs. proper length I for fixed mass M, in units with 2M = 1, for ci = 
0, 0.1, 0.7, and 1.9. In the GR case ci = this is the Einstein-Rosen bridge. For ci = 0.1 the 
radius flares out to infinity so quickly that the code used to make the plot halted at small radius. 
With increasing ci the throat widens, the flare-out inside is slower, and the proper length to the 
curvature singularity increases, becoming infinite for ci > 3/2. 

To more fully compare the Schwarzschild and aether solutions we plot together in Fig. H] 
the radius r{l) and the norm of the Killing vector ^/N{1) for the two solutions with the 
same value of the total mass M. At the internal singularity the norm of the Killing vector 
goes to zero, and the Killing vector is tangent to the constant r surfaces, so the singularity 
is a "would-be" Killing horizon. The surface gravity of the horizon is given by d^/N / dl at 
the horizon. The behavior is easily found using the previous formulae: 

dVN/dl = {l/2){N/By/\Y/r) oc (F - y_)2{i+i^+)/(2+n). (52) 

The exponent is positive, so the derivative of the norm vanishes a.s Y ^ Y^. The surface 
gravity is therefore zero, so the would-be Killing horizon is extremal. 
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FIG. 4: Plot of r(Z) and the norm of the Killing vector \/N{l) for GR and for ci = 0.5, for the 
solution with the same value of the total mass M, in units with 2M = 1. In GR N vanishes at the 
bifurcation sphere at the center of the Einstein-Rosen bridge. In the ae-theory solution the Killing 
vector remains timelike at the throat, but at the internal r = oo curvature singularity both the 
norm and its slope vanish, indicating the presence of a singular extremal Killing horizon. 

3. Charged dust interpretation 

We argued at the beginning of this section that the C2 and C3 terms in the action do not 
contribute to the field equations in the case of static, spherical symmetry and static aether, 
and the C4 term can be absorbed into a simultaneous shift of ci and C3. This enabled us 
to reduce the general case to the one with only ci non-zero. We can nevertheless choose to 
include a C3 term in the action, and by so doing re-express the content of the field equations 
in an interesting way. In particular, if we choose C3 = — Ci, then the Ci and C3 terms combine 
to make (ci/2)F(jbF"^, where Fab — ^aUb — '^bUa- This is just the Maxwell Lagrangian for a 
vector potential Ua, up to a constant factor. We have been treating the contravariant vector 
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M° as the independent field variable, but in this Maxwell-like case it is natural to adopt 
instead the covariant vector Ua as independent. This change just amounts to an invertible 
field redefinition, hence yields the same equations of motion when the metric and aether field 
equations are taken together. With this choice of field variable the theory with C3 = — Ci and 
C2 = C4 = looks quite similar to Maxwell theory, the only difference being the constraint 
term \{g"''^UaUh — 1) in the Lagrangian. It was shown in [a] that this is equivalent to the 
Einstein-Maxwell-charged dust system, with a constant charge to mass ratio fixed by ci, and 
restricted to the sector in which there exists a gauge choice for which the vector potential 
is parallel to the dust 4- velocity. (This is a real restriction.) 

It must be possible to interpret our strange static aether wormhole solution as a charged 
dust solution, but it is at first hard to see what could stabilize the dust unless it is extremally 
charged, which corresponds to the case ci = 2. In fact it was argued in Ref. Q], invoking 



prior results [25| for the charged dust problem, that there is no static solution with static 
aether and general ci. However, our result here shows that argument cannot be correct. The 
contradiction is resolved by the observation that the prior results invoked in Ref. Gj apply 
only if the dust mass density is positive, whereas in the solution at hand this mass density 
turns out to be negative. (In Ref. [26!] the charged dust system with negative mass density 
and constant charge to mass ratio was studied (among other cases) in static axisymmetry, 
and it was shown that every harmonic function determines a solution. Presumably among 
these solutions is the static aether solution found here.) This is related to the negative 
energy density that we already inferred above must be present. With a negative mass 
density, the dust is gravitationally repulsive, and though it has the same sign charge it is 
electrically attractive. (In Newtonian terms, a force F produces an acceleration F/m which 
is opposite to F when m is negative.) Thus the gravitational and electric forces exchange 
their usual roles. The fact that the dust does not then just collapse on itself is perhaps due 
to the associated cost in (positive) electric field energy when the dust is squeezed together. 
Although static, the solution is not regular, since there is an internal singularity. 

4- Solutions with Y G (y_,y+) 

So far we have discussed the solution for Y in all ranges except for Y G (YL, Y^). In this 
range the metric function B is negative, so dr is timelike. Therefore, in order for the metric 
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signature to be Lorentzian, dt must be spacelike. But in this static solution the aether is 
parallel to dt (cf. f|T6l) ). so cannot be timelike. Hence there is no Lorentzian solution with 
timelike aether corresponding to this range of Y . 

As a mathematical curiosity, if we allow > 0, so there are two timelike dimensions, 
there would still be a further restriction for a real solution, since the ratio in is negative 
and raised to the power — / (2 + y+). In order for to be real and positive this power must 
be an even integer m, which implies v = (1/4)(1 — 1/m^). The ratio in fl35l) is also negative, 
and is raised to the different power 1/(2 + which can not also be an integer since is 
not an integer. However, the integration constant rmin can be complex, thus balancing the 
phase of the right hand side of (15^ . and yielding a real solution with signature (+H ). 



V. STARS 



In this section we investigate the static solutions to the theory that, unlike the singular 
wormhole have a regular origin. We first show that there are no such regular pure aether 
solutions, and then turn to the case of fluid stars. 



A. Nonexistence of pure aether stars 

Spherically symmetric self-gravitating "solitons" appear in a number of field systems 
coupled to gravity. For example there are boson star solutions in the Einstein-Klein-Gordon 



system |27l | and Einstein- Yang-Mills theory possesses the Bartnik-McKinnon solutions [28|, 
29|. It is therefore natural to ask whether such "aether star" solutions might exist in the 
vacuum Einstein- Aether theory. The static aether solution studied in the previous section is 
the unique solution with static aether (remember, this means aether aligned with the timelike 
Killing vector), and does not have a regular origin. Thus the only way a regular aether star 
might exist is if the aether has a radial component. We now examine this possibility and 
show that it cannot occur. 

The analysis of Section IIIII showed that local solutions around a general r = tq are 
characterized by three free parameters which may be taken to be A'{rQ), a(ro) and a'(ro). 
If we apply this result at the origin tq = the parameter freedom is restricted. Spherical 
symmetry implies that at the origin the radial component of the aether vanishes, 6(0) = 0. 
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The normalization constraint fll2l) therefore fixes a(0) = e""^''/^, and the r-derivative of 
this constraint imphes that A'{0) and a'{0) are not independent, but rather are related by 
A'{0) = — 2a'(0)/a(0). Thus there is a one parameter family of solutions regular at the 
origin. 

These solutions cannot be asymptotically fiat, for the following reason. The asymp- 
totically fiat boundary condition discussed in Section [1111 would require fixing the one free 
parameter, leaving a unique solution. However, pure ae-theory is scale invariant, so there 
must be at least a one parameter family of solutions much like the Schwarzschild solutions 
of different mass in GR. (By contrast, Einstein- Yang-Mills theory is not scale invariant, and 
admits a discrete family of soliton and black hole solutions.) We conclude that no regular 
aether stars exist. This conclusion was verified empirically by integrating out from the origin 
with different initial parameters, and attempting unsuccessfully to tune to an asymptotically 
flat solution. 

B. Fluid Stars 

Although there are no regular vacuum aether stars, globally regular solutions exist in the 
presence of a static, spherically symmetric perfect fluid with no aether couplings. The fluid 
stress tensor appearing in the metric fleld equation ([5]) is 

Tal = (p(r) + Pir))v,v, - Pir)g,, (53) 

where = e~^^'^{dt)^ is the fluid 4-velocity, p(r) it's mass density, and P(r) it's pressure. 
The metric fleld equation and the Bianchi identity together imply that the sum of the 
aether and fluid energy-momentum tensors is divergenceless. In addition, since the aether 
does not couple directly to the fluid, its stress tensor is independently divergenceless when 
its fleld equation and unit constraint are satisfled. Therefore the fluid stress tensor is also 
independently divergenceless in any solution. Thus, an appropriate system of equations for 
the aether plus fluid case is the (i) metric fleld equation, (ii) aether fleld equation, (iii) radial 
component of V"T^^ = 0, which is the hydrostatic equilibrium equation for the fluid 

P' + lA'{p + P) = 0, (54) 

and (iv) an equation of state p = p{P)- 
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The count of parameters that can be freely specified when integrating out from a regular 
origin is the same as for the vacuum aether solution, except for the addition of a central value 
for the pressure Pq = P(0). For a fixed central pressure, there is just one parameter which 
can be tuned to obtain an asymptotically flat solution. Now there is no contradiction with 
scale invariance, since the central pressure sets the scale and determines the total mass of the 
solution. That these asymptotically fiat star solutions have a static aether can be inferred as 
follows. The field equations can be integrated out from the origin with the aether assumed 
static. The pressure drops to zero at some value r = R, the surface of the star, where the 
static interior solution can be matched to the static vacuum aether solution discussed in 
section HVl This solution is asymptotically fiat, so it must be the unique asymptotically fiat 
solution whose existence is indicated by the parameter count. 

A possible worry about the preceding argument is whether the matching to a vacuum 
solution at the surface of the star could be non-unique. Provided a(r), A{r), B{r) and their 
derivatives are continuous aX r = R one can simply continue the integration into the vacuum 
region. We have not studied this behavior in detail, but it seems that as long as this timelike 
surface is not a characteristic surface of the ODE's corresponding to a spherical wavefront 
of the spin-0 modes of the theory, the field equations will imply the continuity condition. 

To determine a star solution one can fix a central pressure and numerically integrate the 
tt and rr metric field equations and the hydrostatic equilibrium equation flMI) . using the 
equation of state, from the origin to the radius R where the pressure vanishes. There A'{r) 
is continuous so one can use it to match to the vacuum solution. The total mass M can be 
read off from ^ together with ([351), using the definition (ET]), = RA'{R). The area 

of the 2-spheres in such a star solution is strictly increasing as r increases from zero to the 
surface of the star where P vanishes. At that point P' < 0, so according to flMl) A' > 
(assuming positive fluid energy density p). Thus (pTj) implies F > 0, which means that we 
always match to the static aether solution outside of the minimal area 2-sphere. A "throat" 
never occurs in such a star solution. 

1. Constant density stars 

To get a sense of the nature of the static aether star solutions we consider here the simplest 
example, stars with constant energy density interior. Although this does not closely describe 
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FIG. 5: Total mass vs. central pressure Pq in ^ constant density star, in units with p = 1 and 
SirG = 1, for several values of ci. The GR curve asymptotically approaches M = (47r/3)(24/9)^/^ . 
As ci grows the maximum mass decreases, and the curve develops a sharp local maximum and a 
shallow local minimum. 

realistic stars, it is adequate for indicating the behavior of maximum mass limits and the 
stability properties of equilibrium configurations. The discontinuity in the mass density at 
the surface entails via the field equations a jump in A", but A' remains continuous so can 
be used to match to the vacuum solution as described above. 

Graphs of M versus the central pressure Pq for the equilibrium configurations are dis- 
played in Figure O In GR the mass asymptotes to a maximum value as the central pressure 
goes to infinity. Physically, an infinite central pressure would be required to maintain equi- 
librium for a greater mass. As ci increases in the Einstein- Aether case the maximum mass 
limit decreases, and the mass curve develops a local maximum and a very shallow local 
minimum that is only apparent for larger values of Ci. For sufficiently large Ci a second local 
maximum occurs. (We have not attempted to determine the behavior at arbitrarily high 
pressures and for ci approaching 2. Perhaps the series of maxima and minima continues.) 
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FIG. 6: Total mass vs. R for a constant density star, in units with p = 1 and 8irG = 1, for Pq up to 
300 and several values of ci. For small pressures all the curves increase uniformly. For ci = 0.005 
the slope is nearly the same as in GR. For ci = 0.05 a maximum occurs for a large central pressure. 
By ci = 0.6, the first maximum occurs at much smaller pressures, and there is also a minimum. 
For ci = 1.6 a second maximum has appeared. 

The presence of stationary points in the mass versus pressure curves is an indication 
that the stabiUty character of the equiUbrium may be changing {scj. The connection with 
hnearized stabihty arises as follows. If the squared frequency of a mode is positive then the 
corresponding perturbation of the star is oscillatory, while if it is negative the perturbation 
grows exponentially in time. In the borderline case of zero frequency the mode has zero 
energy, hence corresponds to a variation between two static solutions with the same mass. 
Transitions between stability and instability therefore occur at extrema of the mass versus 
R plot shown in Figure [6], where a small displacement of R does not change the mass to first 
order. 

In the GR limit there is no critical point. The mass increases monotonically with central 
pressure, as seen in Fig. so there is no onset of instability. In ae-theory even fiat space is 
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not necessarily stable. The conditions on the q for which all linearized plane wave modes 
have positive squared frequency were found in [IGj]. For example in the pure ci case they are 
—2 < Ci < 1. If we assume the values of Cj are such that very small mass stars are stable, 
then instability can only set in at a critical point of the mass function. As Ci grows larger, 
the curves in Fig. O exhibit extremal points corresponding to the local maxima and minima 
of Fig. [5l At the maximum mass the lowest mode becomes unstable. At the following local 
minimum another zero frequency mode occurs, corresponding to the next mode becoming 
unstable. (It cannot be the lowest mode becoming stable again, since R is increasing with 
increa sing central pressure, implying the presence of a node in the corresponding radial 
mode 30].) Therefore, beyond the maxima in shown in Figure O constant mass density stars 
in the Einstein- Aether theory are unstable. For small values of Ci the central pressure has 
to be very large compared to the density for the star to reach the instability, implying a 
violation of the dominant energy condition. For Ci = 1 the pressure at the onset of instability 
is about 1.28 times the density. 



VI. DISCUSSION 



In this paper we analyzed the Einstein- Aether theory assuming stationary spherical sym- 
metry. We determined the number of free parameters in the corresponding solutions to the 
field equations, and classified the asymptotically fiat ones. The vacuum solution with static 
aether, i.e. aether aligned with the timelike Killing vector, was found analytically up to 
inversion of a transcendental equation. It has a wormhole-like structure, with a minimal 
2-sphere and a singular internal area-infinity which lies at finite affine parameter along a 
radial null geodesic and finite proper distance if Ci < 3/2. The static fluid star solutions 
were classified and studied numerically for the case of constant density, and the maximum 
mass and stability properties were determined. 

Several directions for further work are suggested by these results. Foremost from an 
observational point of view might be to determine the maximum mass limits for a realistic 
neutron star equation of state. For astrophysical applications it would also be necessary to 
determine the structure of rotating solutions. Another open issue is the condition on the 
Ci required by stability of small mass stars. A related question is whether the pure static 
aether wormhole solution is stable. 
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The solutions studied here may be of some help on the unresolved question of energy pos- 



itivity in ae-theory. 
solutions is positive 



.t is known for what ranges of the coefficients q the energy of linearized 
l3, and it is known that for the Maxwell-like special case (03 = — Ci, 
C2 = C4 = 0) nonsingular negative energy initial data exist 12|, |3l|]. The positive mass static 
aether solutions on the other hand have positive energy despite having everywhere negative 
energy density and an interior singularity, as explained in section IIVB2I It would be in- 
teresting to see whether a positive energy result can be established for nonsingular static, 
spherically symmetric star solutions. 

Finally, we found that there is a two parameter family of static, asymptotically flat 
vacuum solutions. Stellar exteriors for different masses form a one-parameter family. The 
other parameter pertains to the radial tipping of the aether away from the Killing vector. In 
a black hole solution the aether must tip, since it cannot remain timelike and be aligned with 



the null Killing vector on the horizon. In the companion to this paper 22| we determine the 
structure of black hole solutions in Einstein- Aether theory. 

Acknowledgements 

This research was supported in part by the NSF under grant PIIY-0300710. 



[1] D. Mattingly, "Modern tests of Lorentz invariance," Living Rev. Rel. 8, 5 (2005) 
|arXiv:gr"^/050209 7 1 . 

[2] CM. Will and K. Nordtvedt, Jr., "Conservation Laws and Preferred Frames in Relativistic 
Gravity. I. Preferred-Frame Theories and an Extended PPN Formalism," Astrophys. J. 177, 
757 (1972); K. Nordtvedt, Jr. and CM. Will, "Conservation Laws and Preferred Frames 
in Relativistic Gravity. II. Experimental Evidence to Rule Out Preferred- Frame Theories of 
Gravity," Astrophys. J. 177, 775 (1972); R.W. Hellings and K. Nordtvedt, Jr., "Vector-metric 
theory of gravity" Phys. Rev. D7, 3593 (1973). 

[3] See, for example, M. Gasperini, "Classical Repulsive Gravity And Broken Lorentz Symmetry," 
Phys. Rev. D 34, 2260 (1986); "Singularity prevention and broken Lorentz symmetry", Class. 
Quantum Grav. 4, 485 (1987); "Repulsive gravity in the very early Universe", Gen. Rel. 



25 



Grav. 30, 1703 (1998); and references therein. 

[4] V. A. Kostelecky and S. Samuel, "Gravitational Phenomenology In Higher Dimensional The- 
ories And Strings," Phys. Rev. D 40, 1886 (1989). 

[5] M. A. Clayton and J. W. Moffat, "Scalar- Tensor Gravity Theory For Dynamical Light Veloc- 



ity," Phys. Lett. B 477, 269 (2000) | arXiv:gr-qc/9910112| . 
[6] T. Jacobson and D. Mattingly, "Gravity with a dynamical preferred frame," Phys. Rev. D 64, 

024028 (2001) |arXiv:gr-qc/0007031 1. 
[7] N. Arkani-Hamed, H. C. Cheng, M. A. Luty and S. Mukohyama, "Ghost condensation and a 

consistent infrared modification of gravity," JHEP 0405, 074 (2004) [arXiv:hep-th/0312099| . 
[8] B. M. Gripaios, "Modified gravity via spontaneous symmetry breaking," JHEP 0410, 069 

(2004) [arXiv:hep-th/0408127| . 
[9] R. Bluhm and V. A. Kostelecky, "Spontaneous Lorentz violation, Nambu-Goldstone modes, 

and gravity," Phys. Rev. D 71, 065008 (2005) fa rXiv:hep-th70412320 1 . 
[10] C. Heinicke, P. Baekler and F. W. Hehl, "Einstein-aether theory, violation of Lorentz invari- 

ance, and metric-affine gravity," Phys. Rev. D 72, 025012 (2005) |arXiv:gr-qc /0504005| . 



[11] V. A. Rubakov, "Phantom without UV pathology," |arXiv:hep-th/0604153 



[12] C. Fling, T. Jacobson and D. Mattingly, "Finstein-aether theory," in Deserfest, eds. J. Liu, 



M. J. Duff, K. Stelle, and R. P. Woodard (World Scientific, 2006) |arXiv:gr-qc/0410001 
[13] C. Fling and T. Jacobson, "Static post-Newtonian equivalence of GR and gravity with a 

dynamical preferred frame," Phys. Rev. D 69, 064005 (2004) |arXiv:gr-qc/0310044] . 
[14] M. L. Graesser, A. Jenkins and M. B. Wise, "Spontaneous Lorentz violation and the long-range 



gravitational preferred-frame effect," Phys. Lett. B 613, 5 (2005) [arXiv:hep-th/0 501223]. 
[15] B. Z. Foster and T. Jacobson, "Post-Newtonian parameters and constraints on Einstein-aether 

theory," Phys. Rev. D 73, 064015 (2006) [arXiv:gr-qc/0509083| . 
[16] T. Jacobson and D. Mattingly, "Einstein-Aether waves," Phys. Rev. D 70, 024003 (2004) 

|arXiv:gr-qc/04 02005| . 

[17] F. A. Lim, "Can We See Lorentz- Violating Vector Fields in the CMB?," Phys. Rev. D 71, 

063504 (2005) fa rXiv:astro-ph/0407437] . 
[18] J. W. Elliott, G. D. Moore and H. Stoica, "Constraining the new aether: Gravitational 



Cherenkov radiation," JHEP 0508, 066 (2005) [arXiv:hep-ph/050521ll . 
[19] C. Fling, "Energy in the Einstein-aether theory," Phys. Rev. D 73, 084026 (2006) 



26 



|arXiv:gr-qc/050 7059]. 

[20] S. M. Carroll and E. A. Lim, "Lorentz-violating vector fields slow the universe down," Phys. 



Rev. D 70, 123525 (2004) fa rXiv:hep-th/040714 9 1 . 
[21] B. Z. Foster, "Radiation damping in Einstein-aether theory," Phys. Rev. D 73, 104012 (2006) 



|arXiv:gr-qc/060 2004| . 

[22] C. Eling and T. Jacobson, "Black holes in Einstein-aether theory," Class. Quantum Grav., to 

appear, (arXiv:gr-qc/0604088 1 . 
[23] B. Z. Foster, "Metric redefinitions in Einstein-aether theory," Phys. Rev. D 72, 044017 (2005) 



[arXiv:gr-qc/0502066|. 



[24] R. M. Wald, General Relativity, (University of Chicago Press, 1984). 

[25] U. K. De and A. K. Raychaudhuri, "Static Distribution of Charged Dust in General Relativ- 
ity," Proc. Royal Soc. A 303 47 (1968). 

[26] W. B. Bonnor, "Equilibrium of charged dust in general relativity," Gen. Rel. Grav. 12 453 
(1980). 

[27] P. Jetzer, "Boson stars," Phys. Rept. 220, 163 (1992). 

[28] R. Bartnik and J. Mckinnon, "Particle - Like Solutions Of The Einstein Yang-Mills Equations," 

Phys. Rev. Lett. 61, 141 (1988). 
[29] M. S. Volkov and D. V. Gal'tsov, "Gravitating non-Abelian solitons and black holes with 

Yang-Mills fields," Phys. Rept. 319, 1 (1999) [ar Xiv:hep-th/9810070] . 
[30] J. B. Hartle, Gravity: An Introduction to Einstein's General Relativity, (Addison Wesley, 

2002), Chapter 22. 

[31] M. A. Clayton, "Causality, shocks and instabilities in vector field models of Lorentz symmetry 



breaking," arXiv:gr-qc/0104103 



27 



